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Abstract. The wave- vector (q) and doping dependences of the magnetic energy, iron moment, 
and effective exchange interactions in LaFeAsO, BaFe2As2, and SrFe2As2 are studied by 
self-consistent LSDA calculations for co-planar spin spirals. For the undoped compounds, 
the calculated total energy, -B(q), reaches its minimum at q corresponding to stripe anti- 
ferromagnetic (AF) order. In LaFeAsO, this minimum becomes flat already at low levels of 
electron-doping and shifts to an incommensurate q at 8=0.2, where 8 is the number of additional 
electrons (8 > 0) or holes (8 < 0) per Fe. In BaFe2As2 and SrFe2As2, stripe order remains stable 
for hole doping down to 8 = —0.3. Under electron doping, on the other hand, the -E(q) minimum 
shifts to incommensurate q already at 5=0.1. 



Iron pnictides have attracted great interest of both experimentalist and theoreticians after 
the discovery of superconductivity with T c =27 K in F-doped LaFeAsOi_ x F x pQ and with even 
higher T c of 38 K in oxygen- free, potassium-doped Bai_ ?; Kj / Fe2As2[2j. Both families of iron 
pnictides have a quasi two-dimensional (2D) tetragonal crystal structure, in which FeAs layers 
are separated by either LaO or Ba layers. The Fe ions form a square lattice sandwiched between 
two As sheets shifted so that each Fe is surrounded by a slightly squeezed As tetrahedron. At 
about 150 K, both stoichiometric parent compounds undergo a structural transition at which 
the symmetry of the lattice lowers to orthorhombic [3J0]. In BaFe2As2 and SrFe2As2, stripe 
antiferromagnetic (AF) order of 0.4-0.8 /xeFe moments [H[S]> aligned ferromagnetically (FM) 
along the shorter 6-axis and antiferromagnetically along the longer o-axis and the c-axis [SJ E] , 
sets in at the same temperature as the structural transition. Transition to a phase with the same 
stripe AF order occurs also in LaFeAsO, but at a 20 K lower temperature than the structural 
transition [7J. Electron doping of the FeAs layers in LaFeAsOi_ x F x suppresses the structural 
and magnetic transitions in favor of superconductivity already at x=0.03 [8j. Also hole doping 
in Bai_j / K y Fe2As2 suppresses the structural and magnetic transition [2j but at much higher 
doping of ~0.15 holes per Fe [9], whereas AF fluctuations are observed up to as high K content 
as y ~0.5. Electron doping caused by Co substitution in Ba(Fei_ rr Co x )2As2, on the other hand, 
suppresses the magnetic transition at x ~0.06. 

Since in many of the superconducting iron pnictides the highest T c is observed at the doping 
level at which the magnetic transition is suppressed, the superconductivity seems to be closely 
related to magnetism and understanding the magnetic interactions between Fe moments in these 
compounds is of utmost importance. Although the theoretically calculated Fe magnetic moment 
and the stabilization energies of different magnetic solutions depend strongly on the employed 
computational method and exchange-correlation functional all band structure calculations that 



the stripe AF order is the magnetic ground state in both parent compounds \10 \ lilt ^ \ 113j. 

Results of spin-spiral calculations and an analysis of effective exchange interactions between 
Fe moments in electron doped LaFeAsO and hole doped BaFe2As2 were presented in our previous 
paper [14] . In this paper we show that electron doping of FeAs layers in M(Fei_ x Co x )2As2 
(M=Ba, Sr) destabilizes stripe AF order as efficiently as in LaFeAsO. 

Self-consistent calculations for co-planar spin spirals in LaFeAsO and MFe2As2, with M=Ba 
or Sr, were performed within the LSDA using the linear muffin-tin orbital (LMTO) method 
in the atomic sphere approximation [15!. The effect of doping was simulated by using the 
virtual crystal approximation, with the doping level 5 defined as the deviation of the number of 
valence electrons from corresponding values for undoped compounds normalized to the number 
of Fe atoms in the unit cells. Positive and negative values of 5 correspond to electron and hole 
doping, respectively. Calculations for LaFeAsO were performed for 5=0.1, 0.2, and 0.3, which 
corresponds to F content of x=S in LaFeAsOi_ x F x . Potassium doping of y=0.2, 0.4, 0.6 in 
Bai_ :y K 2/ Fe2As2 and Sri_j / K< / Fe2As2 was modeled by <5=— 0.1, —0.2, —0.3 (y = 2\5\). Finally, 
electron doping of FeAs layers in M(Fei_ 3; Co ;E )2As2 caused by Co substitution was studied for 
5=0.1. 

The calculations for all doping levels were carried out for experimental room temperature 
crystal structures of the undoped compounds. The tetragonal (PA/nmm) unit cell with 
a=4.0353 A, c=8.7409 A, z La =0.14154, and z As =0.6512 was used for LaFeAsO pQ. MFe 2 As 2 
band structures were calculated for the body centered tetragonal (I4/mmm) unit cell with 
a=3.9625 A, c=13.0168 A, and z As =0.3545 for M=Ba @] and a=3.9243 A, c=12.3644 A, and 
2 As =0.36 for M=Sr [IB]. 

As explained in details in [T3], LMTO calculations for LaFeAsO place Fe d xy bands too 
close to d^^-derived bands as compared to results obtained with the full potential linear 
augmented plane wave (LAPW) method [17] . Better agreement between two band structures 
can be obtained by adding in LMTO calculations an on-site shift of —150 meV to the Fe d xy 
states. In the following, the results calculated for LaFeAsO with the shifted Fe d xy states are 
shown. 
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Figure 1. The q dependencies of the Fe 
magnetic moment (a) and the total energy 
(b) calculated for LaFeAsO with 5=0 (red 
• ), 0.1 (magenta T), and 0.2 (green ■). 



Spin spiral calculations were performed on the base of the generalized Bloch theorem (18j 



under the assumption that the direction of the magnetization is constant within each atomic 
sphere. In these calculations the angle 9 between Fe magnetic moments and the c axis is set to 
90 degree, i.e., all Fe moments lie in the ab plane. The angle <j) between the Fe moment and the 
a axis depends on the wave vector q of a spiral and is given by <pi = + q • R, where R is the 
translation vector. The initial phases for two Fe sites at the positions tj in the tetragonal 
unit cell are fixed by 4> = q • t j . 

The total energy of spin-spirals was calculated for q varying along T (0,0) - X (ir/d,0) - M 
(ir/d,ir/d) - T lines in the two dimensional magnetic Brillouin zone (BZ), where d is the distance 
between Fe nearest-neighbors in the ab plane. At X and M points collinear AF stripe- and 
checkerboard-like spin structures are formed. The calculations were performed for FM (q z =0) 
and AF (q z = n/d±) alignments of Fe moments in adjacent Fe layers, were q z is the z component 
of the wave vector and d± is the interlayer distance. 




Figure 2. The q dependencies of the Fe magnetic moment (a, c) and the total energy (b, d) 
calculated for BaFe2As2 (a, b) and SrFe2As2 (c, d) with AF alignment of Fe moments along the 
c axis (q z = w/d±) and 5= — 0.1 (blue A), (red • ), and 0.1 (magenta ▼). The corresponding 
curves calculated for undoped compounds (6=0) with q z =0 are also shown (red O ). 

By comparing the total energies calculated for q z =0 and q z = ir/d± we found that the 
interlayer magnetic coupling in LaFeAsO is extremely weak. The q dependencies of the Fe 
magnetic moment and the total energy per Fe, E(q), calculated for LaFeAsO with q z =0 and 
5=0, 0.1, and 0.2 are shown in figure [TJ The minimum of E(q) for the undoped compound is 
found at the X point, i.e., for stripe AF order. This, together with the appearance of a local 
minimum along the M-f line, suggests that effective interactions between the nearest (ji) and 
next-nearest (j'2) Fe neighbours are antiferromagnetic with j'2 > ji/2 |14j . The energy of spin 
spirals rapidly increases with the decrease of |q| and and at short wave vectors the self-consistent 
solution becomes nonmagnetic. As a FeAs layer in LaFeAsO is doped with electrons (5 > 0) to 
simulate F doping in LaFeAsOi-^F-r, another local minimum develops at an incommensurate 



wave vector q m i n along the X— M line. E(q m i n ) becomes lower than E(X), i.e., the collinear 
stripe-like ordered AF solution becomes unstable, at 5 >0.1. 

In contrast to LaFeAsO, the energy of spin spirals calculated for MFe2As2 with AF alignment 
of Fe moments along the c direction is significantly lower than for FM one. The corresponding 
curves for the undoped compounds are compared in figure El The change of magnetic order in 
Fe chains running along c from AF to FM costs at the X point 3.5 meV/Fe in BaFe2As2 and 
8.6 meV/Fe in SrFe2As2- The q dependencies of the total energy for both undoped MFe2As2 
compounds are similar to E(q) for LaFeAsO, with the minimum at X and a local one along 
the M-f line. Although the stabilization energy of the stripe-ordered solution, \E(X) — E(0)\, 
is somewhat lower in SrFe2As2 than in BaFe2As2, the energy difference E(M) — E(X) of 38 
meV/Fe in the former is larger than in the latter (29 meV/Fe). This indicates that also j± and 
j2 couplings within a FeAs layer are stronger in SrFe2As2 than in BaFe2As2- 

As it was shown in |14j . when BaFe2As2 is doped with holes, stripe AF order remains stable up 
to the doping level of 5 ~ —0.3, which corresponds to K content of about 0.6 in Bai_ 2/ K J/ Fe2As2. 
The stabilization energy of the stripe-ordered solution is, however, strongly reduced by hole 
doping. The same tendency is also obtained from spin spiral calculations for SrFe2As2 with 5 <0, 
which is illustrated by 5 = —0.1 curves in figure (2) Calculations performed for electron-doped 
MFe2As2 show that the total energy minimum shifts from X to incommensurate q already for 
5=0.1. These results agree with the experimental phase diagram of Ba(Fei_ x Co3;)2As2 according 
to which the structural and magnetic transition at T ~140 K is suppressed at x « 0.06 [19]. 

Recently, strongly anisotropic nearest neighbour couplings along AF (ji a ) and FM (jib) Fe 
chains have been obtained from linear response calculations [20] . We estimated ji a , jib, and j'2 
from a list-square fit to -E'(q) curves calculated using the magnetic force theorem. Preliminary 
results show that the anisotropy of ji a and jib is much weaker than reported in Ref. [20] . 
However, the anisotropy strongly increases if the fit is performed in a small region of q-space 
around the X point. The results of these calculations will be published elsewhere. 
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